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We use a semiclassical approximation to investigate density variations and dipole oscillations of 
an interacting three-component normal Fermi gas in a harmonic trap. We consider both attractive 
and repulsive interactions between different pairs of fermions and study the effect of population 
imbalance on densities. We find that the density profiles significantly deviate from those of non- 
interacting profiles and extremely sensitive to interactions and population imbalance. Unlike for a 
two-component Fermi system, we find density imbalance even for balanced populations. For some 
range of parameters, one component completely repels from the trap center giving rise a donut 
shape density profile. Further, we find that the in-phase dipole oscillation frequency is consistent 
with Kohn's theorem and other two dipole mode frequencies are strongly effected by the interactions 
and the number of atoms in the harmonic trap. 



I. INTRODUCTION 

The recent experimental progress of trapping and cool- 
ing of atomic gases leads to an opportunity for a detail 
study of exciting many body physics pQ. The key to 
this exciting opportunity is that the high experimental 
controllability of these systems. For a dilute atomic mix- 
ture, the quantum degenerate regime can be achieved by 
cooling the system to ultra-cold temperatures. In these 
ultra-cold atomic systems, one of the fascinating control 
parameters has been the two-body scattering length of 
atoms. Tuning the scattering length, simply by applying 
a homogenous magnetic field (Feshbach resonance 0), 
the interactions can be controlled dramatically. The fun- 
damental differences between fermions and bosons take 
place in the quantum degenerate regime. For the case 
of two-component Fermi systems, if the temperature is 
low enough, weakly bound molecules could be formed 
at positive scattering lengths and these molecules can 
undergo Bose-Einstein condensation. By sweeping the 
magnetic field, the scattering length can be further in- 
creased to a divergency and can be changed in sign. In 
this strongly interacting limit, the bosonic molecules con- 
tinuously transform into Bardeen, Cooper and Schrief- 
fer (BCS) superfluid pairs. The physics of this so-called 
BCS-BEC crossover region is relevant to the systems like 
high Tq superconductors, superfluid 3 He, and neutron 
stars. The recent realization of three-component degen- 
erate Fermi gases has opened a new research frontier in 
ultra-cold atoms. These three-component systems have 
relevance to high density quark matter, neutron stars, 
and universal few body physics. Therefore, superfluid 
phases and phase separation in a three-component Fermi 
gas can be treated as an analogous to color superconduc- 
tivity and baryon formation in quantum chromo dynam- 
ics. 

In a recent experiment by Bartenstein et al [3], us- 
ing radio frequency spectroscopic data, molecular bind- 
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FIG. 1: Schematic atomic cloud shapes of a population imbal- 
anced three-component Fermi gas in a spherically symmetric 
harmonic trap. Depending on the interaction between differ- 
ent pairs of fermions and the population imbalance, one com- 
ponent can completely repel from the trap center as shown in 
panel (a). The other two components form a spherical shape 
as shown in panels (b) and (c). See FIG. [4] for quantitative 
details. For a population balance mixture, inner radius of 
the donut shape cloud is zero and we find density imbalance 
through out the harmonic trap. 

ing energies, scattering lengths and the Feshbach reso- 
nance positions of the lowest three channels of 6 Li atoms 
have been determined. As there are three broad s-wave 
Feshbach resonances for the three lowest hyperfine states 
of Li system, one can prepare the system at various 
atomic interactions between different pairs. In another 
recent experimental development by Ottenstein et al [4], 
a three- component Fermi gas has been cooled down to 
the quantum degenerate regime. In this experiment, 
the three lowest hyperfine states of 6 Li atoms are pre- 
pared at a temperature T = 0.37Tp, where Tp is the 
Fermi temperature. During the experiment, the scatter- 
ing lengths for the respective three Feshbach channels are 
varied by sweeping a magnetic field and collisional sta- 
bility of the lowest three channels of 6 Li atoms has been 
studied. The most recent experiment [5] reports the first 
measurements of three-body loss coefficients for a three- 
component Fermi gas at temperature T = 0.59Tp. By 
further lowering the temperature, paring of three com- 
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ponents can be activated. This promising experimen- 
tal observations of three-component superfluidity is yet 
to be realized in near future. Motivated by the realiza- 
tion of these normal degenerate three-component Fermi 
systems [3H5], we study the ground state normal den- 
sity profiles and dipole density oscillations of a three- 
component Fermi gas trapped in a harmonic potential at 
zero temperature. We use a semiclassical theory known 
as Thomas Fermi functional approach to investigate the 
normal state density profiles and dipole density oscilla- 
tions. In general, at high enough temperatures, the parti- 
cles are classical so that the semiclassical approximation 
is valid. In the present paper, we are considering a large 
number of Fermi atoms. These atoms fill every energy 
levels up to the Fermi energy Ep and Ep is larger than 
the ground state energy. Therefore, a three-component 
normal Fermi gas can be well described in the semiclas- 
sical approximation even at very low temperatures pro- 
vided that the number of atoms present in the system is 
large. 

For an equal population system, we find that the den- 
sity imbalanced cloud form a three-shell structure. De- 
pending on the interactions between pairs, density of the 
one component shows non -monotonic behavior as the 
cloud extends from center to the edge of the harmonic 
trap. For an unequal population system, we find very 
similar behavior as the population balanced system, how- 
ever, one component completely repels from the trap cen- 
ter. This repulsion depends on both interaction and the 
population imbalance. The generic shapes of the atomic 
clouds in a harmonic trapping potential for this case are 
schematically shown in FIG. [T] In the presence of har- 
monic trapping potential, we derive expressions for in- 
phase and and out-of-phase dipole oscillation frequen- 
cies as functions of interactions and atom numbers in 
different spin states, we find that the in-phase oscilla- 
tion frequency is independent of the interactions and the 
atom numbers while the out-of-phase oscillation frequen- 
cies strongly depend on them. 

The paper is organized as follows. In the following 
section, we present the basic formalism required for the 
densities and dipole oscillation frequencies. In section 
III, we derive the equations for Thomas Fermi density 
profiles and stability condition and then present our self 
consistent solutions based phase diagram and the density 
profiles. In section IV, we present dipole oscillation fre- 
quencies as a function of interactions and atom numbers 
followed by the discussions and summary in section V. 



II. BASIC FORMALISM 



spectively. Here Ti is the ith cartesian component of the 
position vector r. Let us denote the density of the a 
component is n a and the interaction strength between a 
and v components is U av . 

As we are considering a dilute atomic system at zero 
temperature, and the interactions are short-range in na- 
ture, s-wave scattering channel is dominated over the 
other scattering channels. Therefore, we can neglect the 
higher wave scatterings and the interaction can be spec- 
ified by a single parameter a CT „, which is the s-wave scat- 
tering length between component a and v. The func- 
tional form of the interaction is U av — AnTi 2 a al/ /m. In 
order to derive the ground state densities of the normal 
state at zero temperature, we write the grand canon- 
ical energy as a functional of densities, E = E g — 
J2a=i HvNo- Here the number of atoms in each spin 
state N a = j d 3 rn a and the total energy E g — J ed 3 f, 
where e is given by 
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and A = (3/5)(6tt 2 ) 2 / 3 . Notice we introduced three La- 
grange multipliers (chemical potentials) to constrain 
the number of atoms in each a state. In section III pre- 
sented below, we derive Thomas-Fermi equations for the 
densities from the variation of the total energy. 

In the presence of harmonic trapping potential, the 
lowest lying collective excitations are the dipole or the 
center of mass oscillations. By generalizing the scaling 
method used in Ref. [5], we derive expressions for the 
dipole oscillation frequencies for a three- component sys- 
tem. We assume that the atomic cloud is displaced lin- 
early along r 3 direction. The time varying density of 
the a component is parameterized by the collective co- 
ordinate r) a (t) as n a (r,t) — (rj , r 2 , r 3 — where 
n^\r) is the ground state equilibrium density distribu- 
tion of the er component. Substituting n a (r, t) into total 
energy expression, the variation of the total energy up to 
the quadratic order in r\ is given by 
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We assume that the three-component Fermi gas is 
trapped in a harmonic potential given by V(r) — (ojj_rf + 
u±r 2 + UJ zf 3 )/(2m), where uj z , and m are transverse 
and axial trapping frequencies and mass of an atom re- 
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The dipole oscillation frequencies (u>d) are determined by The classical equation of motion in the matrix form is 
the eigenvalues of the classical equations of motion for rj. given as, 
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III. THOMAS-FERMI EQUATIONS FOR 
DENSITIES 



Before we derive the Thomas-Fermi equations, we in- 
troduce three dimensionless variables; p a = n^^ 3 , r 2 = 

yjh/mui is 



f.2 , ~2 



i 3 and a av = (4tt) a av /C, where £ 
the effective oscillator length with ui = (or^o^) 1 / 3 . The 
dimensionless cartesian component of the vector f are 
defined as fi = [2Trm/ (Tibj)u) 2 ] 1 / 2 ri. In terms of dimen- 
sionless variables, now the trapping potential has a sym- 
metric form given by V[r) = fiuif 2 /(47r). In terms of 
dimensionless parameters, the number equation and the 
total energy are given by 



N a =4n f 2 p G dr, 
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where the dimensionless chemical potential = 
47r p a I (huj) . The Thomas-Fermi equations for the densi- 
ties p a are derived from the variation of the total energy; 
dE/dp a = 0. 
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For the spatial density variations, these equations must 
be solved with the stability condition. The stability con- 
dition can be derived from the determinant of second or- 
der variation of the energy functional, \d 2 Ejdp a dp v \ > 0, 
where d 2 E/dp a dp l , is a 3 x 3 matrix. In terms of dimen- 
sionless variables, this stability condition reads, 
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Similar set of equations for an repulsively interacting 
two-component normal gas in three dimensions are de- 
rived in Refs. [7H9]- For a two-component normal gas in 
two dimensions, the Thomas Fermi equations are derived 
in Ref . [TO] • We have investigated the stability condition 
given in Eq. ([8]). Unlike two-component Fermi gases, 
three-component Fermi mixtures are stable for large den- 
sities and interactions [J]. For a dilute system, we do not 
expect three body collisions to take place so that the 
system will be stable for a wide range of interactions. 

Before we investigate the density variations, let us con- 
struct the phase diagram for a homogenous system. Wc 
consider a set of representative parameters by fixing the 
dimensionless chemical potentials of components one and 
two as pi — 190 and p 2 = 178 (these values are corre- 
sponding to the central chemical potentials of the system 
with particle numbers in the order of 10 5 ). Further, we 
chose the interactions to be 023 = 2.5 and di 3 = 3.6. 
Then by solving the equations Q and ([8]), the phase di- 
agram is constructed in 0,12 - (p 3 —p 2 ) parameter space as 
shown in FIG. [2j When there is no interaction between 
component one and two {a\ 2 — 0), the phase boundaries 
can be easily evaluated analytically. Let us define the 
possible phases by the notation i = I, II, and III, where 
i is the phase with component i. The mixed phases are 
defined as sum of i's (For example, I+II is the phase 
where both component one and two exist). The phase 
boundary between I+II and I+II+III phases is given by 

3/2 3/2 

P-3 = 023P 2 + P13P1 , where /3 CT! , = 3^^/(10^). The 
boundary between I+II+III and II+III phases is given by 

h = (Mi//Si 3 ) 2/3 +VM2 - P1P23/ /3i 3 (P2P23-hl3 23 / Pis)- 
The phase boundary between II + III and III phases is 
given by p 3 = (p, 2 / f3 23 ) 2 ^ 3 . Even though, we have re- 
stricted ourselves to the repulsive interactions between 
two-three and one-three pairs, one can consider attrac- 
tive interactions between these pairs and construct the 
phase diagram in similar manner. 

For density variations, we simultaneously solve equa- 
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FIG. 2: Phase diagram for fixed Hi = (4-7r/?kj)190, /j.2 = 
(47r/ftw)178, fi 2 3 = 2.5, and 013 = 3.6. The labels I, II, and 
III are the phases with component one, two, and three respec- 
tively. Label U denotes the unstable region. 



tions ([5]), ([7]), and ^ for given atom numbers in each 
spin states and given interactions between different pairs, 
fn figure [3j we plot the density profiles for various inter- 
action strengths for an atomic mixture with equal num- 
bers in each three spin states. As seen in the figure, all 
three components exist in the trap center. However, the 
density imbalance remains finite through out the trap. 
The most pronouncing feature is the non-monotonic den- 
sity variations of one component when both repulsive 
and attractive interactions are present. As a result of 
the spatially inhomogeneous trap, the balanced mix- 
ture of atomic cloud forms a 3-shell structure similar to 
the case of population imbalanced two-component gases. 
While the inner most shell contains all three components, 
the outer most shell contains only a single component. 
Even though, we consider a population balanced mix- 
ture in FIG. [3j we find density imbalanced throughout 
the trap, similar to a case of population imbalanced two- 
component system. It should be noted that the density 
imbalance in three- component system here is not due to 
the population imbalance but due to the different inter- 
actions between different pairs. 

The density profiles of an atomic mixture with unequal 
populations are shown in FIG. [4] We find very similar 
density profiles as the case of equal population, however, 
depending on the number of atoms in each spin states and 
the interaction between different pairs, one component 
can completely repel from the trap center giving rise a 
donut shape density profile. This donut shape density 
profile is an unique feature of a three-component normal 
gas in the presence of population imbalance where as this 
feature is absent in a two-component gas. 

For the data in FIG [3] and FIG. [IJ we choose inter- 
action parameters unsystematically to show the different 
behavior of density profiles. Even though the density pro- 
files are shown as a function of scaled parameter r, the 
generic behavior in real space must be the same. This 



can be easily understood if we consider the isotropic case 
where lj±_ = uj z = uj. For this case, the atomic cloud is 
isotropic and the density variations can be presented as 
a function of r, just like in FIG [3] and FIG. |4j 

In experiments, real-space density distributions of each 
components can be obtained by using in situ imaging. 
After trapping and cooling, the system can be prepared 
at any desirable population. As experimental time scales 
are smaller than the spin relaxation time of the atoms, 
fixed spin population can be maintained throughout the 
experiment. Then by sweeping an uniform magnetic 
field, scattering lengths between different pairs can be set 
at different values. As some three-component atomic sys- 
tems (for example, lowest hyperfine spin states of 6 Li g]) 
possess three broad s-wave Feshbach resonances, by tun- 
ing scattering lengths, the interactions between different 
pairs can be varied. However, in a current experimental 
setups, interaction strengths, (i.e, the scattering lengths 
between different pairs) cannot be controlled indepen- 
dently whereas the number of atoms in each spin states 
can be varied independently. 

Density profiles of an repulsively interacting two- 
component atomic gas in a harmonic trap are discussed 
in Ref. [7J [HI EH [H] in the context of collective ferro- 
magnetism. While the authors in these references dis- 
cuss spontaneous magnetization and phase separations 
of two hyperfine spin components, Duine et al dis- 
cuss the nature of the ferromagnetic phase transition in 
the mean field description and beyond mean filed the- 
ory of a homogenous system. The study in Ref. [5] is 
similar to the study in Refs. J7J [T2] , however, authors in 
Ref. [S] go beyond mean field theory and study the effect 
of trap anisotropy on magnetization. Except Ref. |12) . 
all above studies are restricted to two-component atomic 
gases. Salasnich at el [H] investigate the density profiles 
of a three-component gas, however, the authors have re- 
stricted themselves to an equal and repulsive interaction 
between different pairs of fermions. As we have relaxed 
their constraints in the present study allowing to have 
both repulsive and attractive, as well as different inter- 
actions, we find qualitatively and quantitatively very dif- 
ferent results. 



IV. DIPOLE OSCILLATIONS 

By solving for the eigenvalues of Eq. [1J we find the 
three dipole oscillation frequencies, w = uj z and uj± = 



tJi — C ± vB, where 



C= \ Y, K ™( N ° + N ») 



(9) 
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FIG. 3: Density profiles of the three Fermi components for various interactions between pairs. We fixed the number of atoms 
in each spin state to be N a ~ 4.5 x 10 5 . The dimensionless scattering lengths between pairs are: (a) a\ 2 = 1, 013 = — 1, and 
S23 = 4, (b) 2i2 = 0.1, <zi3 = — 1, and 023 = 4, (c) 012 = —2, £13 = — 1, and 623 = 4, (d) 012 = 1, 013 = 2, and 023 = 2. Black 
line: component 1, dashed gray line: component 2, and dotted line: component 3. 
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FIG. 4: Density profiles of the three Fermi components for various interactions between pairs and different numbers. The number 
of atoms and the dimensionless scattering lengths between pairs are: (a) 5i2 = 1, 5i3 = 2, and 023 = 6.6, N\ ~ 3.0 x 10 5 , 
N 2 ~ 3.0 x 10 5 , and N 3 ~ 1.6 x 10 s (b) a 12 = -1, a 13 = 3.6, and a 23 = 2.5, iVi ~ 6.7 x 10 s , N 2 ~ 6.0 x 10 s , and iV 3 ~ 1.9 x 10 5 
(c) 012 = 0.462, 013 = 1.04, and a 2S = 4.51, Ni ~ 2.8 x 10 5 , N 2 ~ 4.7 x 10 s , and N3 ~ 6.7 x 10 s (d) 012 = -0.1, 013 = -0.2, 
and a23 = —0.3 Ni ~ 2.7 x 10 s , ^2 — 1.4 x 10 5 , and N 3 ~ 8.8 x 10 4 . Black line: component 1, dashed gray line: component 2, 
and dotted line: component 3. 
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X - ]T K^K^NvN^ - N a N M - N a N v - (10) 

The parameter K aL , = A av / {rnN a N v ). The eigenvalue 
ujq corresponding to the eigenvector r\\ = r/2 = 773 and 
represents the in-phase oscillation frequency of the sys- 
tem. This mode follows the Kohn's theorem and it 
is independent of interactions and the atomic popula- 
tion. The other two eigenvectors are r)\/Nx = i]2/N 2 = 
— 773/JV3 and r)i/N\ = —^2/^2 = ^3/^3- These two 
modes are corresponding to the oscillations where two 
components are in in-phase while the other one is in out- 
of-phase oscillations. 

In terms of dimensionless parameters, K av is given by 

N a N„ J df df y ' 

In order to see the qualitative behavior of the out-of- 
phase oscillations, we restrict ourself to the equal pop- 
ulation and equal interactions. For the case of equal 
population Ni = N2 = N3 = N and equal interactions 
S12 = 0,13 = o-23 = a, the out-of-phase modes are degen- 
erate and given by 

where p^ = p± + p2 + p^ ■ This density oscillation 
frequency as a function of interaction is shown in FIG. [5] 
for the case of N = 6.5 x 10 6 . For this case case, the 
size of the Fermi gas becomes larger as the interaction 
increases, and then the density oscillation frequency Wd 
monotonically decreases. 

As it has been suggested recently [T3] , the out-of phase 
dipole mode can be excited via two-photon Bragg scat- 
tering in experiments. As the two-photon Bragg scatter- 
ing spectrum is related to the dynamic structure factor, 
measuring the density response spectrum with available 
experimental techniques |14j allows one to measure the 
out-of phase dipole mode frequencies discussed here. 

V. DISCUSSIONS AND SUMMARY 

Even though, we have restricted ourselves to the zero 
temperature in this paper, we do not expect qualitative 
changes at finite temperatures. Even with the attractive 
interaction between some pairs, the system will be in 
normal state if the population imbalance is large. In the 
presence of attractive interactions between different pairs 
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FIG. 5: Out-of-phase dipole oscillation frequency as a func- 
tion of interaction (5) for a symmetric system with N = 
6.5 x 10 6 number of atoms in each spin state. We set the 
interaction between different pairs of fermions to be equal. 

of fermions (and at low population imbalances) , it is pos- 
sible to form s-wave superfluid phases which require equal 
densities. Unlike a two-component Fermi gas, a three- 
component Fermi gas can undergo competitive BCS pair- 
ing and provide non-trivial order parameters. In the pres- 
ence of a superfluid phase, superfluid region will spatially 
phase separate into a concentric shell surrounded by den- 
sity imbalance normal shells. Even though, the experi- 
mental observations of trionic bound states and super- 
fluid phases are yet to be observed in near future, there 
are series of theoretical studies on three-component su- 
perfluidity [IMS]. 

Another promising perspective of a three- component 
Fermi system is to observe universal three-body quantum 
physics called, Efimov states. In the present work, we 
have neglected the appearance of an infinite number of 
three-body bound states (Efimov states) and two-body 
pairing between components. As there is no three-body 
interactions, three-body recombination is forbidden in a 
dilute two-component Fermi gas. 

In summary, we use a density functional approach to 
investigate the Thomas Fermi density profiles and dipole 
oscillation frequencies of an interacting three-component 
normal Fermi gas at zero temperature. We constructed 
a phase diagram for a homogenous system at selected 
representative values of parameters to show that the sys- 
tem possesses a rich phase diagram, even neglecting the 
possible pairing between atom pairs. Even for an equal 
population atomic mixture, we find density imbalance 
through out the harmonic trap due to the different in- 
teraction between different pairs. For an unequal pop- 
ulation atomic mixture, we find that one of the com- 
ponent completely repel from the trap center showing 
a donut shape atomic cloud surrounded by two concen- 
tric spherical clouds. In both cases, density profiles show 
dramatic deviations showing non-monotonic density dis- 
tributions as opposed to the distributions seen in a two- 
component gas. Another qualitative difference from the 
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two-component gas is that the density imbalance even for 
a population balanced normal gas. This density imbal- 
ance in the presence of population balance is due to the 
different interactions between different pairs. We also de- 
rive in-phase and out-of-phase dipole oscillation frequen- 
cies to show how interaction and population imbalance 
effect on these mode frequencies. 
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